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Abstract 

Single-electron transistor with ferromagnetic outer electrodes and nonmag- 
netic island is studied theoretically. Nonequilibrium electron spin distribution 
in the island is caused by tunneling current. The dependencies of the mag- 
netoresistance ratio 5 on the bias and gate voltages show the dips which are 
directly related to the induced separation of Fermi levels for electrons with 

different spins. Inside a dip 5 can become negative. 
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Magnetoresistance of tunnel structures is currently an attractive topic for both experi- 
mental and theoretical studies (see, e.g., Refs. J5~§1)- For tunnel junctions made of ferro- 
magnetic films, the difference as high as 26% at 4.2 K (up to 18% at room temperature) 
between the tunnel resistances for parallel and antiparallel film magnetization has been ob- 
served [TJJ, that allows their application as magnetic sensors. The low temperature values 
agree well with the theoretical result [p]] AR/R = 2P\P2/{1 + P1P2) where Pi and P2 are 
spin polarizations of tunneling electrons in two films, that proves the good achieved quality 
of junctions. As an example, the polarization is about 47% for CoFe, 40% for Fe, and 34% 
for Co H. 

With the decrease of the tunnel junction area, the single-electron charging becomes 
important leading to new physical effects. The study of tunnel magnetoresistance in this 
regime is a rapidly growing field []8| pT4|| (see also Refs. and ||). For example, the enhance- 
ment of the magnetoresistance ratio AR/R due to Coulomb blockade has been discussed 
in Refs. PJI^. The magnetic field dependence of the tunneling current between two Co 
electrodes via the layer of few-nm Co clusters has been measured in Ref. JTO . In Ref. |T2 



the model of small ferromagnetic double tunnel junction (i.e. single-electron transistor |T5j 
(SET-transistor) without gate electrode) has been considered, in which the tunnel resis- 
tances of junctions are different for parallel and antiparallel magnetizations of electrodes, 
thus changing the current through the system. A similar SET-transistor has been also stud- 



ied theoretically in Ref. ||14|| . The very interesting effect of magneto-Coulomb oscillations 
in SET-transistor has been observed and explained in Refs. Strong magnetic field H 

causes the Zeeman shift of two spin bands and their repopulation. Since in ferromagnets 
the densities of states of these bands are different, the Fermi level moves with H leading to 
magneto- Coulomb oscillations. 

In the present letter we consider a SET-transistor which has ferromagnetic outer elec- 
trodes and nonmagnetic central island (see inset in Fig. |l]a). When the coercive fields of two 
ferromagnetic electrodes are different, the standard technique of the magnetic field sweep- 
ing (see, e.g., Ref. |2[]) easily allows to obtain parallel or antiparallel polarizations of outer 



electrodes. In the first approximation the current through SET-transistor does not depend 
on these polarizations because the island is nonmagnetic (the Zeeman splitting is negligible 
because of small H). However, if the electron spin relaxation in the island is not too fast (es- 
timates are discussed later), then the tunneling of electrons with preferable spin orientation 
creates the nonequilibrium spin-polarized state of the island (similar to the effect discussed 
in Refs. [|16| , |17|1 in absence of the Coulomb blockade). This in turn affects the tunneling in 
each junction and leads to different currents I p and I a through SET-transistor in the parallel 
and antiparallel configurations. 

We will calculate the dependence of the relative current change 5 — (I p — I a )/I p on 
the bias and gate voltages (we call 5 magnetoresistance ratio despite for finite voltage this 
terminology could be misleading). Nonzero 5 is already the evidence of the nonequilibrium 
spin state in the island. Moreover, the voltage dependence of S shows the dips, the width of 
which directly corresponds to the energy separation between Fermi levels of electrons with 
different spins in the island. 

We consider SET-transistor consisting of two tunnel junctions with capacitances C\ and 
C 2 . Induced background charge Q as usual [|T3J describes the influence of the gate voltage 
(general equivalence relations for finite gate capacitance are discussed, e.g., in Ref. [0). 
We assume that the voltage scale related to the magnetic polarization of ferromagnetic 
electrodes [|] and the voltage scale of the barrier suppression [19| are large in comparison 
with the single-electron charging energy (that is a typical experimental situation). Then 
the polarization of outer electrodes can be taken into account by the difference between 
the tunnel resistances R^ 2 and Rf 2 for electrons with "up" and "down" spins. The total 
junction resistances R± = (1/R% + l/Rf)' 1 and R 2 = (I/-R2 + V-^D -1 do not depend on 
the magnetic polarizations Pi and P 2 of electrodes, while "partial" resistances are given by 
R% = 2R i /(l + and Rf = 2^/(1 - i*) (similar to the model of Ref. [§). Notice that 
Pi describes the polarization of tunneling electrons || which is different (typically even in 
sign) from the total electron polarization at Fermi level (the latter one determines, e.g., the 
period of magneto-Coulomb oscillations ||). 
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We assume that the energy relaxation of electrons in the island is much faster than the 
spin relaxation. So, we characterize the nonequilibrium spin state by the difference AEp 
between Fermi levels for "up" and "down" spins while both distributions are determined 
by the thermostat temperature T. The spin diffusion length is assumed to be much larger 
than the island size (that is a typical experimental situation - see Ref. ]rB|), so the spin 
distribution is uniform along the island. 

The equations of the "orthodox" theory for single-electron transistor |7|Jl5| (we assume 
Ri ^> Rk = h/e 2 ) should be modified in our case. The energy gain for tunneling to 

(+) or from (-) the island through ith junction is different for "up" and "down" electrons, 

W^n) = W? =F WtHn) = W? ± ^ (1) 



s 



T(ne + Q )T(-1 NJ 



a 2. 

Here n is the number of extra electrons on the island (as usual the electron is assumed 
to have positive charge e), Cs = C\ + C2, and V is the bias voltage. The corresponding 
tunneling rates satisfy usual equation [7| 

r ±(n) = ?£W (2) 

#R? (l - eM-WiHn)/T)) ' 

where s = u,d denotes spin. The average current / through SET-transistor can be calculated 

as 

/ = ]>>[n+H-rrH]a(n), (3) 

n,s 

where a{n) is the stationary solution of the master equation [Q 

da(n) jdt = H n ± !) r r(^ ± 1) - °( n ) r f{ n )}- ( 4 ) 
Finally, the Fermi level separation AEp should satisfy selfconsistent equation 

AE F pv/T = Y.[T u l + {n)-Ti+{n) 

a(n), (5) 



-rr(n)+ir(»j 



where r is the electron spin relaxation time for the island, p is the density of states (per 
spin), and v is island's volume. We introduce also the dimensionless spin relaxation time 
a = r/e 2 pv(Ri + P 2 ). 

The signs of polarizations Pi and P2 can be changed using the external magnetic field, 
that interchanges resistances P" and Rf. So the current I p for the parallel magnetization 
(P1P2 > 0) is different from the current I a when one magnetization direction is reversed, 
P2 — > —P2 (the change Pi — > —Pi obviously gives the same result). Figure |l|a shows the 
numerically calculated dependence of the magnetoresistance ratio 5 (solid line) on the bias 
voltage V for the SET-transistor with parameters C\ = C2, Qo = 0, T = 0, |Pi| = | P^ | = 
30%, and a = 0.1. For the upper curve (shifted up for clarity) we assumed P 2 = Pi 
while R2 = 5Pi for the lower curve. The S—V dependence shows the oscillations with the 
same period e/Cj as for the Coulomb staircase. The existence of oscillations is a trivial 
consequence of the charge dynamics in SET-transistor, similar to the effect discussed in 



Refs. 12.14 



More interesting features seen in Fig. |I|a are the triangular-shape dips near the bias 
voltages 

V=[e/2 + ne + (-l) i Q ]/C i , (6) 

at which the derivative of the I-V curve abruptly increases (the Coulomb staircase for I p 
shown in Fig. |l]a by dashed lines is better seen for the lower curve). The edges of a dip 
correspond to the alignment between the Fermi level in an electrode and one of the split 
Fermi levels for electrons with different spins in the island. Hence, the dip width AV is 
directly related to the Fermi level splitting, = AEpCs/eCi. Somewhat similar effect 
(the finite voltage caused by injected spin current) without Coulomb blockade had been 
observed in the "spin transistor" [IBJ fabricated using the thin film geometry with a size 
scale of 50 /mi; the typical signal scale was only about 30 pV. The small size of the SET- 
transistor island leads to the very strong enhancement of the Fermi level separation (see Eq. 
(D) and the corresponding voltage scale. 



The width of the dips in Fig. [l|a increases with voltage because the larger current provides 



larger AE F (the crude estimate is AE F = aI{P\ — P 2 )eRs). Actually in the case |Pi| = |P 2 | 
shown in the figures, the nonequilibrium spin distribution in the island occurs only for the 
antiparallel magnetization of electrodes; for parallel magnetization AEp = because the 
spin currents to/from different electrodes exactly cancel each other. When |P X | ^ \P 2 \ the 
complete cancellation does not occur, and the dip shape is determined by two different values 
of AEp leading to the trapezoid-like shape instead of the triangular one. 

It is interesting that the magnetoresistance ratio S can be even negative within the dip 
range (see Fig. |l|). This can be understood in the following way. The I a -V curve for the 
antiparallel magnetization generally goes below I p -V curve because the Fermi level splitting 
(which is larger in the antiparallel case) decreases the effective voltage for the spin band 
which provides the easier tunneling. However, this also splits the kinks on the I-V curve 
leading to the increase of I a (and decrease of 5) within the splitting range. For sufficiently 
steep kink, I a can become even larger than I p (negative 5). This also explains why the 
dips are more significant for larger tunnel resistance ratio (see Fig. |I]a) when the Coulomb 
staircase is more pronounced. 

Increase of the spin relaxation time r leads to larger AEp and, hence, increases 5 as 
well as widens the dips, that is illustrated in Fig. |l|b (5 = for r = 0). The change of the 
polarization amplitudes |Pi| and |P 2 | leads to similar effects. Crudely, the magnetoresistance 
is determined by the product a\PiP 2 \, while the exact shape of the 5-V curve depends on 
each of these parameters. 

In the limit of large bias voltage the magnetoresistance ratio can be found analytically 
using the following expression for the current: 



However, the formula for 5 is rather lengthy, so we present here only the result for small a, 



IRx/V= l-(a/2)(P 1 -P 2 ) 2 / 




(7) 



5 = 2a|P 1 P 2 |, 



(8) 
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and the expression 5 = 2a|PiP 2 |/[l + 2a (1 — (\Pi \ — I ) 2 /4)] for the case R\ — i?2- 

The finite temperature smears the features of the 5-V dependence (see Fig. 0a), but 
obviously does not change 5 in the large-bias limit. The dips disappear when T becomes 
comparable to AEp while the oscillations disappear at higher temperatures determined by 
the single-electron energy scale e 2 /Cs- 

Notice that two series of dips determined by Eq. (H) coincide in Figs, [l] and ^|a. With 
the change of the background charge Qq by the gate voltage, these two series will move in 
opposite directions. The dips can be also seen on the 5-Qo dependence which is shown in 
Fig. |^b for different bias voltages V. The dip position moves with V. There are two dips 
per period, however, one of them is much less pronounced because of relatively large ratio 
of tunnel resistances. 

To estimate the parameters of a possible experimental realization, let us assume Co-Cu- 
Co SET-transistor. (Notice that Al-Co-Al SET-transistor has been already fabricated j9], 
however, aluminum is not suitable for our purpose because of its superconductivity.) The 
polarization \P\ = 30% used in figures is a conservative value for Co. The spin relaxation 
rate r for nonmagnetic island, which is the most crucial parameter of the effect, depends 
much on the material quality. In Ref. [20| r ~ 10~ 7 s has been reported for very pure Cu 



at T = 1.4 K (the similar value has been found in Ref. [20f for Al, while r ~ 10~ 8 s have 
been reported for Al in Ref. [0). Let us choose r = 10~ 8 s. Then using p = 9 x 10 21 
eV _1 cm -3 for Cu, = 10 5 f2, and the island volume v = 200nm x 50nm x 20nm, we get 
a = 0.35. Hence, the effect of nonequilibrium spin distribution should be rather strong, 
and we could expect the magnetoresistance ratio S up to ~ 10% (5 is significantly enhanced 
near the Coulomb blockade threshold - see Fig. [I]). This large value allows to consider the 
possible applications of such a device. For C*s ~ 3 x 10~ 16 F the dips of the 5-V dependence 
could be observed at temperatures below ~ 0.2 K while the oscillations could be noticeable 
up to T ~ 1 K. 

In our theory we have neglected the Zeeman splitting ±g^BH/2 because the typical 
coercive fields are relatively small, H ~ 10 2 Oe 0. Hence the corresponding energy scale is 



very small, AE ~ 10~ 6 eV ~ K, and the effect can hardly be observed. 

We have discussed dc case only. If ac voltages are applied to the bias and/or gate 
electrodes, the similar formalism can be used. However, in this case the dynamic solution of 
the master equation (^) should be used instead of the stationary solution, and also the left 
side of Eq. (^) should be replaced by [d(AE F )/dt + AE F /r]pv. The nontrivial dependence 
starts when the frequency / of the applied voltage becomes comparable to r -1 . These 
frequencies are within the experimentally achievable range, so such an experiment could be 
used for the direct measurement of the spin relaxation time r. 

In conclusion, we have considered the SET-transistor consisting of ferromagnetic elec- 
trodes and nonmagnetic island. We have predicted that the nonequilibrium spin distribution 
in the island leads to a considerable magnetoresistance which has a specific dependence on 
the bias and gate voltages. In particular, it shows the dips directly related to the Fermi 
level splitting. 
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FIGURES 

FIG. 1. The dependence of the magnetoresistance ratio 5 on the bias voltage V for (a) two 
different ratios R2/R1 and (b) for several values of the dimensionless spin relaxation time a. Inset 
in (a) shows the schematic of the SET-transistor with ferromagnetic (F) outer electrodes and 
nonmagnetic (N) island, while the dashed lines show the I p -V curves (arbitrary units). 

FIG. 2. (a) The S-V dependence for different temperatures T and (b) the dependence of S on 
the background charge Qo for several bias voltages. 
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